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Abstract 

The equilibration of matter and onset of hydrodynamics can be understood in the 
AdS/CFT context as a gravitational collapse process, in which "collision debris" create 
a horizon. In this paper we consider the simplest geometry possible, a flat shell (or 
membrane) falling in the holographic direction toward the horizon. The metric is a 
combination of two well known solutions: thermal AdS above the shell and pure AdS 
below, while motion of the shell is given by the Israel junction condition. Furthermore, 
when the shell motion can be considered slow, we were able to solve for two-point func- 
tions of all boundary stress tensor and found that an observer on the boundary sees a 
very peculiar quasiequilibrium: while the average stress tensor < 7)i,y> contains the 
equilibrium plasma energy and pressure at all times, the spectral densities of the corre- 
lators (related with occupation probabilities of the modes) reveal additional oscillating 
terms absent in equilibrium. This is explained by the "echo" phenomenon, a partial 
return of the field coherence at certain "echo" times. 
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1 Introduction 



Observation at RHIC of collective flows in relativistic heavy ion collisions, well described 
by ideal hydrodynamics [2, 3] have lead to a paradigm shift in the field toward studies of 
strongly coupled plasmas. AdS/CFT correspondence [1] is one of actively pursuit directions, 
both for understanding of strongly coupled gauge theories in general and properties of 
Quark-Gluon Plasma (QGP) in particularly: for recent review see [4]. 

This paper continues the line of work of our two previous papers, [5] and [6], which we 
will call I and II respectively below. Those papers include extensive introduction explaining 
our approach to the problem, which we will not repeat here. Let us just say that in those 
works we dealt with "elementary" collisions - calculating the shape of the falling string 
between two departing charges and its hologram at the boundary - which in the QCD 
language can be related to e+e~ annihilation into a pair of heavy quarks or pp collisions. 
Thus in I and II there was no horizon on the gravity side and no temperature or entropy 
on the matter side. In this paper we address these issues, related with heavy ion collisions 
and equilibration. 

Properties of equilibrium strongly coupled conformal plasma is by now well studied in 
significant detail in the AdS setting, in static thermal AdS-BH metric suggested by Witten. 
Many details about quasinormal modes of this metric and in particularly the correlators 
of stress tensors and their spectral densities are known, we especially recommend [8, 9, 11, 
10]. Recent developments included flowing near-equilibrium state, with slowly /gradually 
deformed horizons and derivation of hydrodynamics up to second order in gradients: we 
will not use those here, see references in a review [4] . 

The most challenging task of the theory now is the understanding of how matter 
manages to equilibrate so rapidly in RHIC collisions, and what exactly such equilibration 
means microscopically. The success of hydrodynamics in describing RHIC elliptic flow 
data seems to suggest the thermalization time of the order of 0.5 /m/c, yet its mechanism 
remains unclear. The quest for its mechanism involves studies of various phenomenology and 
theoretical approaches. We will not attempt to review it and just mention one approach to 
ideas of which we will refer below, the "Glasma model": see Ref.[7] and references therein. 
It is based on classical Yang-Mills equations and starts with the so called Color Glass 
Condensate initial condition. While the coupling is assumed weak in this approach, strong 
coherent fields make its behavior nonperturbative. 

Our approach also attempts to address the transition from Glasma-like initial co- 
herent gauge fields to incoherent near-thermal QGP, but in strong coupling and thus the 
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AdS/CFT setting. Since the vacuum corresponds to extremal black hole solution - pure 
AdS geometry without a horizon - while the thermal field theory is dual to AdS black hole 
with a horizon, the main issue is how deposition of an extra mass - from "debris" of the 

collision - leads to excitation of the extremal to non-extremal black hole and dynamical 
creation of the horizon. In simpler words, we have to follow some kind of a gravitational 
collapse. 

In I and II we already discuss qualitatively how falling debris created in the collisions 
- e.g. large number of falling strings between departing charges - make a falling matter 
shell or membrane. While in I we studied its fall ignoring its own weight, we have to include 
it now, as its near-horizon breaking is entirely due to adding the mass of the shell to total 
gravity. We have argued in those papers that in principle one can address the problem 
by following the motion of two membranes, the shell and that describing horizon (a la 
"membrane paradigm" see book by Thorne and collaborators [16]). 

For now we do not study motion of those two membranes in realistic geometry, for 
obvious reasons starting with the simplest geometry possible. We assume that the matter 
shell (and thus the horizon) is flat - that is independent on our world 3 spatial coordi- 
nates, and moves only in the 5-th holographic direction*^. The setup of this gravitational 
collapse model is described in Sec.2. The equation of motion for the shell is given by Israel 
junction condition[17] which we solve numerically. We find how its trajectory depends on 
the property of the shell; but in all cases the distant observer sees its slow approach to the 
horizon at late time. 

Early works along this path include important papers [12, 13, 14] which we partly fol- 
low. They consider a collapsing shell geometry, but unlike us they use as a probe some bulk 
scalar field and we don't find their boundary conditions at the shell sufficiently convincing, 
and we tried to improve those by using gravitons instead. 

What the observer sees as the shell falls is dual to the evolution of = 4 SYM 
from certain initial ensemble to the final thermal equilibrium. The insight into the problem 
of thermalization is thus obtained by studying various observables - the induced stress 
tensor and its correlation functions on the boundary - while the shell is somewhere in its 
process of falling. The former is given by one-point function and the metric above the shell, 
which in our geometry is time independent AdS black hole metric. Thus the "single point 
observer" who is only able to measure the average density and pressure would be driven 
to the conclusion that the matter is fully equilibrated at all times. More sophisticated 

*'^the equilibration in this setting is not due to spatial gradient as in hydrodynamics 
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"two point observer"- measuring the stress tensor correlation functions - would however be 
able to see the deviations the from the thermal ensemble. We compute those deviations in 
Sec. 3, using various components of the graviton perturbations to probe the gravitational 

background with a shell. 

As a significant technical advance of this work, we show that unique prescription for 
boundary conditions for the gravity waves on the shell follows from the junction condition 
itself. Thus we show how correctly propagate the graviton wave across the shell, relating the 
obvious infalling conditions near the AdS center to what is seen on the boundary. Explicit 
expressions are obtained for two-point function in near equilibrium stages, when the shell 
is close to the horizon. We solve the wave equations both numerically and using the WKB 
approximation, finding good agreement between the two. Possible implications of the results 
are finally summarized and discussed in Sec.4. 



2 Gravitationally collapsing shell in AdS 
2.1 The background metric 

Our setting includes the basic AdS background, described by the metric ds^ = _ 
Its holographic coordinate z is zero at the boundary (UV) and infinity at the AdS center 
(IR). The problem considered is a simple generalization of Israel's original problem, which 
was collapsing spherical shell in asymptotically flat 3d space. And it shares its main feature: 
although the shell is falling with its radial position time depending, the gravity both inside 
and outside it is time independent. Furthermore, inside a sphere there is no influence of the 
shell's existence at all: the famous statement going back to Newton himself. The gravity 
outside only knows the total shell mass. 

It is not difficult to prove that the same is true for fiat shell in the AdS setting as 
well. Starting with a generic form: 

ds^ = -A{z, t)df + B(z, t)dx^ + C(z, t)dz^ (1) 

one can set S = ^ by a coordinate transformation. The metric has to satisfy the vacuum 
Einstein equation 

G^,^. - kg^y = (2) 

with A = 6 both above#^ z < Zjn and below z > Zm the shell's position Zm- We will also 

*^The reader is reminded that the coordinate z is inversely proportional to radial coordinate, thus some- 
what counterintuitive inequalities. 
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refer to those as "outside" and "inside" regions below. 

The tz component tells us that dfC = 0. The tt,zz equations are used to obtain: 

1 



Ciz,t) 



z^{l + kz^) 



A{z,t) = F{t)^^^ (3) 

F{t) can be dropped by a rescaling in the t coordinate. Now we require the metric should 
reduce to the AdS form infinitely far away from the shell, as z — > 0, we can have k = — ^, 
then outside the shell the metric is in form of translationally invariant AdS-black hole(AdS- 
BH). On the other hand, inside at the AdS center z ^ oo we have to set fc = to suppress 
the term. Therefore the inside is just AdS metric. 

So the background metric is the combination of AdS-BH and AdS, with the two 
metrics separated by a shell. We will use the metric in the usual form 

ds"^ = ^{-f{z)dt'^ + dx'^ + dz^/f{z)) (4) 

4 

with / = 1 — ^ (or/ = 1) outside (or inside) the shell position Zm- x and z are both 
continuous in order for ^ to match. 

Note that a singularity z = Zh is outside the region where the former metric is 
used, as < z^. This does not mean that there is no horizon in the problem: in spite of 
pure AdS metric inside, the dynamical horizon and trapped surface (both time dependent) 
do exists. 

2.2 Israel junction conditions and the falling shell 

As in [5] , the strings in AdS bulk were modeled by a shell (membrane) , the action of which 
is given by: 

Sm = -P J d^a^J-detgij (5) 

where gij is the induced metric on the shell, p is the only parameter characterizing the 
shell. 

We use Lanczos equation to study the falling of the shell: 



[Kij]-9iAK] = -4Si3 (6) 
where [K,,] = K± - Kr^ and K^, = n^i^ + 



5 



We parametrize the induced metric on the shell as follows: 



dsl = -^ + ^ (7) 

We choose a = t,z,x and i = t,x Assume the EOM is given by z(r), t(r). 

The norm is determined from the condition nadx" = and = 1 as: 

na = {-^,0,^) (8) 

Note here the norm points to the AdS center(2; = oo). Therefore we have +:inside; — :outside 
The curvature K is calculated as follows: 

_ Uf' + 2fz _l 

K:cx = -- (10) 

z z 

Sij is determined from the shell action: Sij = -^^^j- = PQij- (6) becomes: 

[Krr-grrK] = '^ (11) 

2 

[Ka^x-9xxK] = -'^ (12) 



(11)^ VrTi^-v/7Ti^=^ (13) 
(14) is solved by (13) with integration constant determined! 



(¥)' + (t7^)'(1-/)'-^ (15) 



6 ' ' ^2kIp' ^' ■" 2 
The falling velocity seen by distant observer is given by: 



dt i 4p A. 3 /-i _ f\ 

6 ^ 2;^y^ J) 



(16) 



Suppose the shell starts falling at z = zq > with vanishing initial velocity. The 
horizon radius z^ can be expressed in terms of zq and n\p: 
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4 = l-/(.o)=4(l-f^)^ (17) 

DO 

Note (17) implies another constraint n^p < 6. The independent parameters Zq and p 
should be estimated from the initial conditions (e.g. energy density, particle number), these 
will determine the equilibrium temperature of the evolution. 

With chosen parameters, it is easy to integrate (16) to give the trajectory of the 
shell. Wc have plotted the shell trajectory in Fig.l. It shows three stages of falling, initial 
acceleration (z = zq + #t^), intermediate near-constant velocity fall, and the final near- 
horizon freezing with exponentially small deviation (z = z^ — e~**). 

Finally, what are the physical meaning of the parameters of our model, the initial 
position zq and the shell tension p? The wave functions of the colliding nuclei are believed 
to be [18] concentrated at the certain momentum scale called the "saturation scale" Qs, 
which depends on collision energy and nuclei: it is about 1.5 GeV at RHIC. The holographic 
coordinate has the meaning of the inverse momentum scale in the renormalization group 
sense, thus its initial value should be identified with inverses saturation scale zo = l/Qs- 
Further "falling" corresponds to motion of the wave function into the infrared direction, till 
equilibrium is reached. The initial temperatures at RHIC (and thus Zh = I/ttT are believed 
to be about .35 GeV, with ttT 1.1 GeV: thus the expected inequality Zh > zq is satisfied. 

(Physical fireballs not only equilibrate but also expand, with T decreasing by a factor 
3-4 in RHIC collisions. This would correspond to a departing horizon toward larger z, as 
suggested in [19], but this feature is of course not included in the present model.) 

The corresponding tension of the shell p may be calculated from (17). Although we 
do not follow this direction here, one may attempt to calculate it in a particular model of 
the collisions. In particular, the so called Lund model prescribe how many color strings 
per transverse area is created, and as the shell is but an approximation to all those strings 
falling together, its tension may be identified with the sum of the tensions of all the strings. 

3 Correlation functions of the Gauge Theory 

With the complete gravity background at hand, we are ready to study the property of the 
dual gauge theory under thcrmalization. We already commented that since the metric is 
asymptoticly AdS-BH, the one-point function of the stress tensor is the same as thermal 
SYM result. Therefore we consider the two-point functions as the expected place to find 
deviations from the equilibrated thermal ensemble. 
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Figure 1: The shell trajectory as a function of time. It starts at rest at z = zq with a 
constant acceleration, followed by a constant falling and eventually approaches the horizon 
in a exponential fashion. The parameter we choose are K^p = 1 and Zh = l.Q 

The most standard way is to study a probe field in the background which is dual to 
some boundary operator, then use AdS/CFT prescription to find the correlation function. 
The simplest probe field used in early works was the bulk scalar [12, 13, 14]. However we 
choose to use various components of the metric field, /im„(m, n = t,x), which is dual to the 
boundary stress tensor. One obvious reason is to avoid the introduction of new field into 
the problem. The other reason involves the matching condition of hmn on the shell, as will 
become clear later. Also we only probe the geometry after the creation of the shell. 

Thus we solve for the metric perturbation hmnit, x) propagating in the bulk specified 
above. This is a very difficult task in general, because the shell is always falling and thus 
there are time dependent boundary condition at the shell: this makes a straightforward 
Fourier decomposition in time impossible. However, a possible simplification can be made 
if the Fourier mode is much faster than the falling of shell, i.e. the condition a; S> ^ holds, 
we may consider the shell as quasi-static. In this limit, we can trust the Fourier mode and 
study the problem in the usual way. In other words, we only trust quantities obtained for 
uj ^ With this argument, we can in principle study two-point function in any stage 
provided the mode is fast enough. 



3.1 Matching Condition on the Shell 

Israel junction condition in general should be applicable for any gravity fields. Therefore, 
one can apply it for background plus graviton perturbation and from the latter obtain the 
matching condition for the graviton. 

Similar as in Sec.2.2, we start with the Lanczos equation: 

[Kij - QijK] = -K^Sij = -nlPQij (18) 

which we cast into a different form: 

[K^,] = ^g., (19) 

The zcroth order (in /imn) Lanczos equation have already been used above, for calculation 
of the trajectory of the shell. Now we require vanishing of the first order terms in Lanczos 
equation 

'f« = "«(a|^ + a^^) (20) 

with a = t,z,x and i = t,x. = (~f jOj f) remains unchanged, so the variation of K^j 
comes entirely from Christoffel: 

5K.,=n^5T%^^ (21) 

We choose the gauge h^z = and further assume hmn = hmn{t,w, z), where x = 
xi,y = X2,w = xs- This will not affect the two-point function because the gravity back- 
ground is rotationally invariant in Calculating the variation of Christoffel to the first 
order in hmn and noting z = 0,i = = (quasi-static limit), we find the only 

non-vanishing components of 5Kij is: 



5Kxy = 


2 1 1'xy 


SKrx = 


— -^dzhtx 


S^xw ~ 


2 '^Z'''XW 


6Krr = 




6Ktw = 


——dzhtw 


SKxx = 


^^dh 

2 ^Z"'XX 


()K — 


^^dh 
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We have omitted some components involving index y: those can be obtained by the substi- 
tution X — >■ y from those hsted above. Plugging to (19), we have: 



"•xy — "■xy 

dzhtx — dzh{x = — 3^^*^ 

— hf i — ^ 

iHx — 'Hx'' ~ ^ 



„ xw 

3z 



hxw — h^w 



dzhtt - '^^Mt = 



Ozlltw — Ozlltw — 3^ 

I, _ 'Hw 

hw-^ 



dzhxx - y/ldzhix = --^Kx 

hxx — h^x 

9zhww ^/~f9zh{i]yj = n~3zhyjyj 

hww ~ ^WW (^^) 

J- 

We use from here on hmn and hmn for metric perturbations inside and outside the 
shell respectively. All the quantities are evaluated on the shell z = Zm- The other identities 
follow from the continuity of induced metric across the shell. 

Note the jump in time coordinate, we have to do the Fourier transform in a consistent 
way*^: / dtout^^'^"'^* ^ VJ ^ dtine^'^/^^^'^ , The indices "in" and "out" represent quantities 
inside and outside the shell. We use oj as the frequency measured by clock outside the shell, 
the corresponding frequency inside is given by Thus we obtain from (22): 



*^The convention we use is hmn{t,w) = J /imn(w, g)e'"* ^'^^dtdw 
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9zHy = jidzhxy + -^h^y) 
Hx = htx 

^zHx = A^idzhx + ^^htx) 



V7 3z 



^/7 



9zHw = jidzKw + 
Ht = ^Jlhtt 
dzMt = dzhtt + -^^tt 

2 



xw I 



y =—h 

"■XX ^j^i^xx 



9zHx = jidzhxx + '^^hxx) 

y - —h 

"■WW ~ 

^zh{)jyj = ~^{dzhyjyj H ^UIU> ) (^3) 

All the quantities are evaluated on the shell z = Zm 

Prom here on, we define u=4 = z'^UTf in accordance with the literature [201. The 
axial gauge is just hmu = 0- The metric perturbations can be classified into three channels, 
according to [9]: scalar channel: hxy shear channel: htx and hxw sound channel include hu, 

htw 5 ^ww and h(nx — hxx ~1~ ^yy 

The metric perturbations satisfy linearized Einstein equation, and they are deter- 
mined up to residual gauge transformation hmn — '^m^n — '^n^m where should 
preserve the axial gauge chosen above. Instead of fixing the gauge, one can look for gauge in- 
variant combination in each channel. The behavior of these gauge invariant objects encodes 
complete information of retarded correlator [9, 21]."^^ 



**there is a subtlety in the sound channel, which will be elaborated later 
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3.1.1 The scalar channel 

The gauge invariant object is simply ^3 = hly outside the shell and ^3 = hxy inside. The 
EOM of 03 is given by(with / = 1 corresponding to ^3) 

The prime denote derivative with respect to u. The matching condition between ^3 and 
can be easily obtained from (23): 



1 

7?^ 



4' = ^i'p3 + ^<p3) (25) 

Besides the matching condition (25) , we also need boundary condition at AdS center 
to uniquely fix the solution of (t)^{(f>3) in the bulk up to normalization. The boundary 
condition we use at AdS center 2; = 00 is infalling wave or regular wave. With these 
conditions, we are ready to proceed. Let us start from inside the shell: ^3 is given by the 
solution to the following equation: 

where fm = f{um)- It origins from the jump of frequency across the shell. The EOM 
can be solved in terms of cylindrical function: 



V Jm 

H^^\2^J^/J-^^) ^<-|g| (27) 
K2(2V9'-a;2//^V^) \q\ > 



VJm 



While outside the shell, (pi^ can be written as a linear combination of two independent 
solutions, which we select to be the infalling wave and the outfalling waves at the horizon*^ 

4 = c+cf>+ + c-cf>^ (28) 

where are solutions to (24). If extrapolated to u > Um, ~ (1 ~ ■u)^*'^/^ as 
u — s- 1. (25) gives: 



*^Note that it is just a formal basis for the solution above, we don't use it below the shell and there is no 
horizon singularity there. 
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c_- 0+p-</,+'gl"="- ^^^^ 

where P = ) (^'3 - ^^3) , Q = ^h- 

We expect to recover the equilibration because, as the shell approaches the "horizon" , 
the region where geometry deviates from AdS-BH shrinks exponentially. The ratio |^ — > 00. 
All deviations from equilibrium should become exponentially small as well, as we expect 
from any other small deviations from equilibrium. 

We would like to confirm this limit by our matching/boundary condition. As the shell 
approaches the "horizon", Um — 1,/m Oj we may disregard the third situation in (27). 
We want to calculate the ratio ^ to the leading order in or (1 — Um)- The correction 
to 03 (tf) = (1 — n)^**^/^ is linear in while the correction from asymptotic expansion of 
Hankel function is of \/Ym. Therefore we may simply use 0^ = (1 — -u)^*"^/^. 

Let's focus on the first situation 03 = H'^\2\y/u), with X = ^ y~ — = -^={1 + 



0{fm))- The correction due to can also be ignored at leading order. Using the asymp- 
totic expansion of Hankel function, we find the leading order result cancels exactly in the 
denominator, while the counterpart survives in the numerator. We end up with: 



The asymptotic ratio (30) tends to infinity universally for any REAL uj, correctly 
recovering the AdS-BH limit. 

3.1.2 The shear and the sound channels 

The gauge invariant object in the former case is (f){ = qhL + uhL outside the shell and 
01 = qhtyj + ^=hxyj inside. Here again the frequency inside is scaled by —k=- The matching 
condition between 0{ and 0i from (23) turns out to be the same as the scalar case, up to 
constant normalization: 



= * 



The EOM of 0{ is given by(with / = 1 corresponding to 
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With enough luck, we note satisfies the same EOM as ^3 (26), which combined 
with (31) guarantees the same asymptotic ratio (30). 

Finally we consider the sound channel. The gauge invariant object is 

(f)l = ^ (^'^q'^Mt + 4w^^L + '^^'^hlw + (g^(2 - /) - io^)hL^ outside the shell and 
02 = 2q%t + -^htw + T^hww + {q^ - j^)haa inside the shell. 

This time wc do not seem to have simple matching condition as (25) and (31). An 
exception is at g = 0, in which case, we have: 



4>2 = 02 

'■>2 = —r^\.92 -^^02) (33) 



For case g 7^ 0, an easy way to avoid general discussion is to take advantage of the 
residue gauge. It can be shown by a proper choice of residue gauge, (33) still holds. The 
particular gauge choice of course does not shift the retarded correlator. We include a brief 
justification for (33) in Appendix.B 

The EOM of 02 is given by (with / = 1 corresponding to 02): 



Jn _ q^f^8q'f + Aq' + 9iV^f-6iv\ f, _ 

uf{q^f + 2q^-3u^) '^^ 
Su'^f - 6u;V - Aq'^u'^uf - 2q^uj'^u + 3a;^n + 2gVu - q^f + q^fu + Aq^f 



MM) 



We again find 02 satisfies the same EOM as 03 (26), which combined with (31) 
guarantees the same asymptotic ratio (30). 

Summarizing the discussion above, we have found the EOM of gauge invariant objects 
for three channels: (24), (32) and (34). We also find universal matching condition for all three 
channels (up to a constant normalization): 



0{ = 0a 
] 



i' = -^{<t>a+^<t^a) (35) 
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where a = 1,2,3 and is given by (27). We may from now on forget about the 
geometry inside the shell and simply consider (35) as a boundary condition at the shell. 

3.2 Retarded Correlators and their spectral densities 

In this section, we will solve for the gauge invariant objects and extract the retarded cor- 
relator. As in [20], we switch from (pa to Za = uchaj {'kTY . which couples to the boundary 
stress tensor. The equation satisfied by Za can be simply derived from (24), (32) and (34): 



Zz —Z2. H 72 — Z3 = (36) 

„ (g;^ - gV)/ - uuj^f , iJ^ - q^f 
Z\ ji-^ Zi -I — — Zi = U (37) 

^„ 3u;Hl + u^) + qH2u^-3u^-3) ^, 
^ w/(3a;2 + g2(^2_3)) ^2 + 

3a;4 + g4(3 _ 4y2 + ^4) + ^2(4^2^2 _ 5^2 _ 4^3 

u/2(3a;2 + g2(^2_3)) ^2 - (38) 

Without a simple analytical expression of the ratio between the infalling and outfalling 
waves, we study numerically the solutions to (36). The boundary conditions from (35) are 
given by: 



where a = 1,2,3. We have written the boundary conditions for the three channels 
in a uniform way. This is achieved by an appropriate scaling in Za and does not affect the 
two-point functions. 

The retarded correlator are obtained according to the prescription specified in [21]: 
Ga = -tt'^N^T^ lim - ha ln(«)^ (40) 



u— >0 \ 2Za 



where ha = — \{q^ — uP')'^. 

The retarded correlators are extracted from the boundary behavior of numerical so- 
lutions to (36). The retarded correlators correspond to some off-equilibrium state. We 
compare them with the counterparts in equilibrium state, which are obtained from numer- 
ical solutions to (36), but with infalling boundary conditions. 
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In particular, we study the retarded correlator of momentum density, energy density 
and transverse stress. They are related to the gauge invariant correlator by [9]: 



1 

2 



Gtx,tx - o ,,2_^2 ^1 (^1) 



Gtt,tt - ^— — 3^G^2 (42) 



Gxy,xy - (43) 



We focus on the spectral density, which is defined by: 



X^l.v,p\{'^, q) = -2ImG,^^^px{uj, q) (44) 

The spectral densities of the transverse stress arc plotted at q = 0, q = 1.5 and 
q = 4.5 in Fig. 2. Each plot includes five values of u^, corresponding to different stage of 
thermalization. The thermal spectral density of the transverse stress Xxy,xy "^Iso included 
for comparison. All the non-thermal spectral densities can be viewed as some oscillation on 
top of their thermal counterpart. The first period of oscillation occurs near oj = g*^. The 
oscillation damps in amplitude and grows in frequency as the medium thermalizes*^, i.e. 
Um approaches 1 from below. This effect is clearly illustrated in Fig.3, where we plot the 
relative deviation R = ^^^'^I'^-^y^-y . 

'Kxy,xy 

Parallel to the case of transverse stress, we also plot the spectral density of momentum 
density, at g = 1.5 and = 4.5 in Fig.4 (The spectral density of momentum density at g = 
vanishes identically). Each plot include five values of ix^, corresponding to different stage 
of thermalization. 

The relative deviation R = jg plotted in Fig.5. We again observe the 

damping of amplitude and growing of frequency as the medium thermalizes. 

The spectral density of energy density, at g = 1.5 and q = 4.5 are plotted in Fig.6. 
(The spectral density of energy density again vanishes). Each plot includes five values of Um, 
corresponding to different stage of thermalization. The relative deviation R = '^"'"th''^"'" 
is plotted in Fig. 7. We find a very sharp peak in the first period of oscillation, which is 
removed from the final plot for a better comparison. We again confirm the non-thermal 
spectral density relaxes to thermal one in the qualitatively the same way as described in 
the previous cases. 



*®Similar behavior is also observed in thermal spectral density[21, 22] 

*'^Here the frequency refers to the oscillation in spectral density. It is understood as the reciprocal of to 
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Fi gure 2: (color online)The spectral density of transverse stress Xxy,xy in 

unit of 7r2iV2r4, at 

q = left, q = 1.5 middle and q = 4.5 right. Plotted are spectral densities at different stages 
of thermalization: black asterisk(nm = 0.1), red box(nm = 0.3), blue circle(nm = 0.5), 
green cross(nm = 0.7), brown diamond(tim = 0.9). The thermal spectral density is also 
included(pink point) for comparison. The parameter we will keep using from here on is 

4p = 1 




Figure 3: (color online)The relative deviation R at q = left, q = 1.5 middle and q = 4.5 
right. Different stages of thermalization are indicated by: black asterisk(?im = 0.1), red 
hox{um = 0.3), blue circle(Mm = 0.5), green cross(iim = 0.7), brown diamond(Mm = 0.9). 
As Uffi approaches 1, i.e. the medium evolves to equilibrium, the oscillation decreases in 
amplitude and increases in frequency, thus the spectral density relaxes to thermal one 
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Figure 4: (color online)The spectral density of momentum density Xtx,tx in unit of n'^N^T^ 
at q = 1.5 left and q = 4.5 right. Plotted are spectral densities at different stages of 
thermalization: black asterisk(um = 0.1), red hox{um = 0.3), blue circle(tim = 0.5), green 
cross(iim = 0.7), brown diamond(nm = 0.9). The thermal spectral density is also in- 
cluded(pink point) for comparison 




Figure 5: (color online) The relative deviation R at q = 1.5 left and g = 4.5 right. Different 
stages of thermalization are indicated by: black asterisk(nm = 0.1), red hox{um = 0.3), blue 
circle(nm = 0.5), green cross(um = 0.7), brown diamond(tim = 0.9). As Um approaches 1, 
i.e. the medium evolves to equilibrium, the oscillation decreases in amplitude and increases 
in frequency, thus the spectral density relaxes to thermal one 
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Figure 6: (color online)The spectral density of energy density Xtt,tt in unit of vr^A^^T^, at q = 
1.5 left and q = 4.5 right. Plotted are spectral densities at different stages of thermalization: 
black asterisk(tim = 0.1), red hox{um = 0.3), blue circle(nm = 0.5), green cross(um = 0.7), 
brown diamond(tim, = 0.9). The thermal spectral density is also included(pink point) for 
comparison 




Figure 7: (color online)The relative deviation R at q = 1.5 left and q = 4.5 right. Different 
stages of thermalization are indicated by: black asterisk(nm = 0.1), red box(nm = 0.3), blue 
circle(nm = 0.5), green cross(um = 0.7), brown diamond(iim = 0.9). As Um approaches 1, 
i.e. the medium evolves to equilibrium, the oscillation decreases in amplitude and increases 
in frequency, thus the spectral density relaxes to thermal one 
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In order to explain the observed phenomenon, we would like to obtain some analytical 
formula for the spectral density. This is possible in the final freezing stage, where there is 
a simple asymptotic ratio (30). To the leading order in 1 — Um, the boundary condition as 
u —>■ 1 is: 



Za = C+Z+ + c^Z- ^ c+(l - n)^ + c_(l - u)^ (45) 

For the purpose of illuminating the problem, it is enough to focus on Z^, the EOM 
of which has the simplest form. Its EOM (36) is solvable in terms of Heun function. 
However the property of Heun function is not fully understood.*^. Wc have to use some 
approximation method, and in the regime of large lu the WKB treatment is appropriate . 
Following [26, 23], we obtain the expression of Z^ near the boundary up to normalization(see 
Appendix. A for details of the treatment): 



(^Hi^\2^uj'^-q-^y/^) - |^e-2^'^«oijf ^(2Va;2 - g^^)) u < -\q\ 



(46) 



+ ^-^^ (\ - |^e-2*'-''o) e-'^'=o 12(2^92 -a;20I) > u; > -\q\ 

where the constants ao, 60) cq are defined by: 



, , /l-.s-(l-n2) 1 , 

hmu^x I duJ _ ^^^^ =ao- -ln{l - u) 







hniu^i I du^j _ ^^^^ = 60 - -ln{l - u) 



uo 



S2(l -U^)-l 



U2)2 



Co 



S = |-| 

to 

no = (47) 

The retarded correlator is calculated from the prescription (40). We have dropped 
an additional contact term: 7r2iV2r''(l/2 - (u;2 - ^2)2^). 



#8, 



see [10] for a discussion 
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{ln{uj — q ) + isgn{uj)iTg{sgn{uj), ao)) s < 1 

4 



^^^^2l!fc^(/n(g2 - cj2) + is5n(w)7re-2|^l«5(s5n(w), 60)) s > 1 



where we have defined g{e, c) = l^tel-^2ZllJ/l'\^ - 

In the large uj limit, the non- logarithmic term in the second case is exponentially 
suppressed. Therefore the lowest order result in c_-/c+ agrees with the zero temperature 
one. [23]*^. The higher order correction is due to the emergence of wave from the "horizon" . 
The phase difference between infalling and outfalling waves gives rise to the oscillating 
behavior in to. We have also calculated (48), with C-/c+ numerically obtained and restored 
the contact term. The result show very good agreement with retarded correlator obtained 
in Sec. 3. 2 in region of large uj. 

The physical interpretation of (48) is most clear at vanishing momentum, s = 0. The 
n-th order correction appears as 

^n^R ^ fMp^2msgn{u:^+^) (^l^e-^-^o^ (49) 
Combined with the ratio of ^ (30), the n-th order correction can be written as 



5^G^ ~ ^4-ngiu)(in(l-«m)-2ao)n ^qq-j 

Recall the WKB solution of the incoming wave at s = 0: 



^ ^ gl-l/2^-iu> f^S'du+iu!t 

liniu^i / S'du = — -/n(l — u) + ao (51) 
Jo 2 

Taking into account the time factor e*'^*, we see — (Zn(l — Um) — 2ao) is just the time for the 
wave to travel back-and-forth in WKB potential. We define the echo time 

techo = i-ln{l - Um) + 2ao) (52) 

in which the wave makes a roundtrip. The n'th order correction to the two-point function 
(spectral density) has an echo time of n * techo, with a suppressed amplitude, obviously the 
n-th echo. 

*^the imaginaxy part has an opposite sign, which is due to a different convention in Fourier transform 
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This resembles the usual echo phenomenon with a sound reflected by a wall. Fur- 
thermore, the oscillations in spectral density as a function of frequency result, after the 
Fourier transform, in a peak as a function of time at the echo time. The peak is a result 
of many harmonics coherently added together: while smooth equilibrium spectral densities 
correspond to (thermally occupied) field harmonics which are completely incoherent to each 
other. We thus interpret "echo" as a partial re-appearance of coherence which was present 
in the original "color glass" fields at the collision moment. As the shell keeps falling to- 
ward the "horizon", Um — > 1, the echo time tends to infinity and the medium loses all the 
coherence. 

Although we derived it in near-thermal position of the shell, the echo phenomenon 
by itself is not restricted to the final freezing stage but exists throughout the whole process: 
its manifestation is in the spectral density of Fig.3, Fig.5 and Fig.7. Looking for "echo" in 
dynamical (time dependent) shell is perhaps worth addressing in later work. 

3.3 Can the bounciry observer see what happens below the shell? 

In the rest of the paper we have chosen the boundary condition to correspond to the infalling 
wave at u-infinity (the AdS center), which leads to the solution (27) and its consequences 
discussed above. This particular choice is natural in the standard setting, when all the 
probes (the source and the sink) sending gravitational wave from the AdS boundary u = 0. 

Now we switch to another issue: the gravitational wave emerging from inside, below 
the shell, u > u„i. The motivation for studying this case is as follows. All stationary black 
hole metrics are such that no signal from the inside the horizon can propagate outside it: 
in particular, the geodesies do not cross horizon. But in the falling shell case we consider, 
the metric coincides with black hole one only above the shell, while inside it is the AdS^ 
solution without the horizon, since gravity of the shell produces no eS^ect inside it. The 
question then is, can an observer on the boundary see what happens inside the shell? 

Thus a wave sent from below would propagate all the way till the shell without any 
problems, and scatter on it. So we are looking now for a solution, which in the region inside 
the shell contains both outfalling and infalling waves, while outside {u < Um) it has only 
outfalling wave propagating toward the boundary. At the shell the matching condition is 
again given by (23), and the precise relation between outfalling and infalling waves can be 
worked out parallel to what we did above. As always, a generic graviton is split into three 
channels, with the same eqns. 

The gravitational wave outside the shell contains only the outfalling component. 
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This is a solution in thermal-AdS metric which, if extrapolated to (non-existing) region 
Ujn < u < 1 would be originating from the horizon. Thus at n = 1 it is ^{ (1 — u)~*'^/^, 
it is understood as the behavior since u < Um- The wave inside the shell can be written as: 

= Cin<i>in + Cout<l>out (53) 

With similar argument as before, we can approximate (j)a = (1 — ti)"*"^/^. Matching solu- 
tions at both sides and according to Israel condition (35), we obtain the asymptotic 
reflection ratio 

Cout 

_ {l/8-4p/6)y'l-ul 

Cout ^ 

which, at the shell approaching the horizon Um 1 tends to zero. This means as the shell 
approaches the horizon, the portion of the wave reflected by the shell disappears, and thus 
all the wave emerging from below the shell is transmitted! 

(To convince the reader that this conclusion is correct, here is another argument. As 
we have shown in the beginning of the paper, as the shell approach horizon one recovers 
the solution for the AdS-BH background (without any shell), with only the infalling waves 
on either side of the horizon without reflection (black horizon). The solution with only 
outfalling waves we are now describing is its complex conjugate.) 

The conclusion that even at late time our collapsing shell is not that black as a 
horizon, as signals from inside the shell can be seen, look surprising and worrisome at first. 
Note however that these signals are both strongly red-shifted and exponentially delayed by 
the warping factor y^JJumj- As — 1, both effects become infinitely strong, in practical 
sense precluding the boundary observer from seeing what happens below the shell. 

4 Conclusion and Outlook 

Continuing the line of research started by our papers I and II, we have built a gravita- 
tional collapse scenario, describing equilibration of conformal strongly coupled plasma in 
the AdS / CFT setting. Using a simplified geometry we approximated falling collision debris 
by a single flat shell or membrane, falling from its initial position (given by the saturation 
scale) to its horizon (given by the equilibrium temperature). The setting itself provides 
inequality between the two scales, satisfied at RHIC. 

The main simplification of the paper is that this shell is flat - independent on our 
world 3 spatial coordinates. Therefore the overall solution of Einstein eqns reduces to two 
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separate regions with well known static AdS-BH and AdS metrics. The falling of the shell 
is time dependent, its equation of motion is determined by the Israel junction condition, 
which we solve and analyzed. We also determined how final temperature (horizon position) 
depends on initial scale and shell tension. 

In statistical mechanics it is known that fluctuations at different scales are indepen- 
dent, and that long-range fluctuations (IR) need more time to be equilibrated than the 
short-range (UV) ones. With the AdS/CFT setting and our geometric simplification, this 
fact is taken to its perfect form. It is connected to the well known gravitational phenomenon: 
gravity of a sphere is independent on its size outside it, and is completely absent inside. 
In the equilibrating gauge theory it means that mean quantities at all scales above some 
z < Zm{t) (corresponding to sliding position of the shell) are exactly as in equilibrium, while 
those below it are absolutely unaffected, being the same as in vacuum without any matter. 

This is quasiequilibrium in the title. More specifically it means the following. In 
this geometry a "single point observer" - who is only able to measure the average density 
and pressure - would be driven to the conclusion that the matter is instantaneously equili- 
brated at all times. However more sophisticated "two point observer" who is able to study 
correlation functions of stress tensors would be able to observe deviations from the thermal 
case. We computed them explicitly, calculating a number of spectral densities at various 
positions of the shell, corresponding (in quasi-static approximation) to different stages of 
equilibration. 

The equilibration process can roughly be divided into three stages: the initial accel- 
eration, intermediate relativistic falling and final near-horizon freezing. By studying the 
graviton probes corresponding to three different combinations of hmn, we calculate the re- 
tarded correlators of T„in on the AdS boundary, The matching condition of hmn on the 
shell is given by a variation of Israel junction condition. In the quasi-static limit, we study 
the retarded correlator of all graviton probes. We have shown that the collapsing geometry 
correctly reproduces the AdS-BH limit: as the shell approaches the "horizon" the ratio 
of the infalling wave and outfalling wave tends to infinity. We further confirmed this by 
numerically study the spectral density for transverse stress, momentum density and energy 
density, which allows us to see deviations between geometry with shell and equilibrated one 
(AdS BH). 

We find that the main deviation between the non-thermal and equilibrium spectral 
densities are some oscillations. As the time goes on and the shell is at the position closer 
and closer to the horizon, these oscillations become exponentially smaller in amplitude 
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and higher in frequency, eventually disappearing in the equilibrium state. In this sense 
we get numerical control over the relaxation process. In the final freezing stage, when the 
membrane is close to the horizon, and in large a; regime, we even find analytical expressions 

for these deviations using the WKB method. 

Oscillations in spectral densities as a function of frequency are further explained by 
the "echo" effect, producing peaks at certain "echo" times in the response functions. We 
expect at those times partial restoration of field coherence which was present at the initial 
time of the collision in system's wave function. For the near equilibrated medium we have 
the echo time analytically, from the WKB solution. The echo time tends to infinity as the 
medium thermalizes. 

The echo phenomenon arises from the phase coherence between infalling and out- 
falling waves in the bulk at certain times: we expect it to exist in all gauge theory with a 
gravity dual. It is also interesting to extend the current study to scalar probe and vector 
probe. One can further study the effect of echo on electromagnetic (e+e~) spectral densities 
related to production spectra [27] which can be observable in collisions. 

An attentive reader will notice that apart of small discussion of upward moving 
waves, we we have not yet addressed the dynamics of the horizon formation, deferred for 
later studies. When this paper were near-completed, we learned about an interesting work 
by Hubeny,Liu and Rangamani [15] who discuss certain null geodesies related observation 
points at the boundary. Since phases of the waves add coherently near geometric optics 
paths (null geodesies), this is another interesting form of coherence, although perhaps un- 
related to our WKB "echos" . 

Finally, one may now think about relaxing our main assumption - flatness of the 
shell, e.g. by including first corrections resulting from slow variations of its position. In 
this case the metric above the shell would become time-dependent, allowing a "single point 
observer" to see some relaxation dynamics as well. 
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A WKB treatment of (36) 

In order to apply WKB method to (36), we need to convert it Schrodinger- type equation. 
Introducing a new field ip = ^ ^^^^3, (36) becomes: 



u;\l-s\l-u^)) -3 + 6n^ _ 
with s = |5l~l-wisa large parameter to justify WKB. There are two singularities 



in (55): u = 0,u = 1. The term proportional to may vanish at uq = if s > 1. We 

discuss two cases separately: Is<l, IIs>l and focus on a; > only. The solution for 
a; < can be obtained easily from the solution for u; > by the substitution c_|- ^ c_ . 
Case I s < 1: Away from the singularities, the WKB solution to (55) is given as: 



with s' = j '-f:^^ 

Near the singularities, (55) becomes: 



n ^ 

< + "^'^-^"^'^ V^ - -^V = ^ V = V^i??^'^'^(2a;^r^V^) 



4(1 — uy 



On the other hand, we have: 



u ^ S' ^ \ - — — f S'du = 2V^Vl - 
\ u Jo 

1 1 

u ^ 1 S' ^ — / S'du = ao- -Inil - u) (58) 

2(1 — u) Jq 2 

Matching the WKB solution with the approximate solutions near the singularities 
(using asymptotic expansion of Hankel function), we obtain: 



26 



l-\-iu) 1 — iuj 

V' = c+(1 — n) 2 +c-[l — u) 2 

~ (^//f^(2u;v'l-s2V^) + '-^e-^''^'"'H^^\2^v^/l - s'^V^^ (59) 

Case.II s > 1: This case is a little more complicated because WKB approximation 
breaks down near u = uq. Away from u = 0, 1, uq, we have the following WKB solutions: 



u>uo V^ = 5'-V2e±*"^«o^'<^" 

u<uo = S'-V2e±-/:°^'<^« (60) 



where S' = ,^^g^,S' - ■ /SiEZEl 



We first match WKB solutions at u = uq- Near u = (55) becomes: 



V'" + a;^a(n - no)V' + btp = (61) 
(61) can be solved by Airy functions: 



Using the asymptotic expansion of Airy functions: (x > 0) 



_2 3/2 

e 3 

Ai{x) ~ 



2V7ra;V4 



2-3/2 

Bi{x) ~ 



sm(|x^/2 _|_ lyr) 



-^Haj) V 1 /d (63) 



V^xV4 

We obtain the following match between WKB solutions: 



27 



f (V^ + ^>) + §(V'^ - - + ^V'^ (64) 

Next we match WKB solutions with approximate solutions near singularities similarly 
as case I. Finally we have: 



tp = C-^{1 — U) 2 +C-[1—U) 2 

/ \ p^CO 



with liniu-ti S'd-u = 6o ~ ^^'^(1 ~ '")>/o"° ^^'dn = cq 
Summarizing all cases, we have: 



(65) 



(66) 



(l - ^e-^-'-o) e:ff!Ii^2(2V?^0I) 
-¥ (l + l^e-^^'^^o) e-'^^/2(2v/^23^V^) < a; < |g| 

(l + l^e-^-'-o) e:if£if2(2V?^V^) 
. (l - ^e-^^'-^o) e-'^«/2(2V9'-w'v^) > a; > 

B Gauge Choice for the Sound Channel 

The aim of the section is to show it is possible render the following matching condition by 
a proper choice of gauge: 



Vf 3« V7 

With the help of (23), (67) can be simplified to: 



4' = ^-h' + i^^ (67) 
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2 — f 

-^hL = 



haa ^^yy ii'aa 

haa' + '^haa - (2 - f)hL' + fh^ = (68) 

where all the quantities are evaluated at n = Um- 

We note the residue gauge degree of freedom implies that it is sufficient to satisfy 
(23) up to a gauge choice. In particular, we could add pure gauge solution to the sound 
channel: /i5«"S"= inside the shell and /i9""se,/ outside. According to [24], the only pure gauge 
solution that touch hia{haa with / = 1) is: 



,y _ V7(l + n^ + 2a;' 



^gauge,, 

** U 
,gauge,f _ -qUJ SXCSUlU - qUJUy/J 



^g^augej ^ ^ arCSm . - (gg^ 



Now, it is enough to satisfy: 



^-^jfUa = 



h'aa + ^ha^ " (2 - f)hU + f'Ha = (70) 

where haa = haa + ^ and hla = h[a + It is easy to see it is always possible to 

satisfy (70) with a proper choice of constants A and B. 
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